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^ '. Abstract 

Simulated dark matter haloes are fitted by self-similar, universal 
II . density profiles, where the scaled parameters depend only on a scaled 

O I (truncation) radius, H = R/vq, which, in turn, is supposed to be inde- 

pendent on the mass and the formation redshift. The further assump- 

03 ■ tion of a lognormal distribution (for a selected mass bin) of the scaled 

radius, or concentration, in agreement with the data from a large 
statistical sample of simulated haloes (Bullock et al. 2001), allows 
(at least to a first extent) a normal or lognormal distribution for other 
scaled parameters, via the same procedure which leads to the propaga- 
tion of the errors. A criterion for the choice of the best fitting density 
profile is proposed, with regard to a set of high-resolution simulations, 
where some averaging procedure on scaled density profiles has been 
performed, in connection with a number of fitting density profiles. To 
this aim, a minimum value of the ratio, \xrf\/crsr^ = 1^ ~ ^*l/'7s»?) is 
required to yield the best fit, where rj is the arithmetic mean over 
the whole set; r]* is its counterpart related to the fitting density pro- 
file; asrj is the standard deviation from the mean; and rj is a selected, 
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scaled i.e. dimensionless parameter. The above criterion is applied 
to a pair of sets each made of a dozen of high-resolution simulations, 
FMOl (Fukushige & Makino 2001) and KLAOl (Klypin et al. 2001), 
in connection with two currently used fitting density profiles, NFW 
(e.g., Navarro et al. 1997) and MOA (e.g., Moore et al. 1999), where 
the dependence of the scaled radius on the mass and the formation 
redshift, may be neglected to a first extent. With regard to FMOl 
and KLAOl samples, the best fits turn out to be MOA and NFW, re- 
spectively. In addition, the above results also hold in dealing with rms 
errors derived via the propagation of the errors, with regard to the dis- 
tributions of scaled parameters. The sensitivity error of simulations is 
also estimated and shown to be less than the related, standard devia- 
tion, that is a necessary condition for detectability of accidental errors. 
Some features of the early evolution of dark matter haloes represented 
by fitting density profiles, are discussed in the limit of the spherical 
top-hat model. Though the related matter distributions appear to 
be poorly representative of simulated haloes, unless the (mean) peak 
height is an increasing function of the mass, the results are shown to 
be consistent, provided considerable acquisition of angular momentum 
takes place during the expansion phase. 
keywords - cosmology: dark matter. 

1 Introduction 

According to a wide number of both analytical and numerical studies (e.g., 
Cole & Lacey 1996; Syer & White 1998; Navarro et al. 1995, 1996, 1997, 
hereafter quoted as NFW97; Moore et al. 1998, 1999, hereafter quoted as 
MOA99; Fukushige & Makino 2001, hereafter quoted as FMOl; Klypin et 
al. 2001, hereafter quoted as KLAOl; Fukushige & Makino 2003, hereafter 
quoted as FM03), dark matter haloes which virialize from hierarchical clus- 
tering show universal density profiles, p = p{r; po,ro), where po is a scaling 
density and tq is a scaling radius. In this view, smaller haloes formed first 
from initial density fluctuations and then merged with each other, or were 
tidally disrupted from previously formed mergers, to become larger haloes. 

The density profile is (i) self-similar, in the sense that it has the same 
expression, independent of time (e.g., FMOl), and (ii) universal, in the sense 
that it has the same expression, independent of halo mass, initial density 



perturbation spectrum, or value of cosmological parameters (e.g., NFW97; 
FMOl; FM03). A satisfactory fit to the results of numerical simulations is 
the family of density profiles (e.g., Hernquist 1990; Zhao 1996): 

^\ro) (r/ro)'^[l + (r/ro)"]x ' ^ a ' ^ ^ 

for a suitable choice of exponents, a, f3, and 7. 

This family includes both cuspy profiles first proposed by Navarro et al. 
(1995, 1996), NFW97, (a,/3,7) = (1,3,1), hereafter quoted as NFW den- 
sity profile, and the so called modified isothermal profile, (a, P, 7) = (2, 2, 0), 
which is the most widely used model for the halo density distribution in anal- 
yses of observed rotation curves. It also includes the perfect elhpsoid (e.g., 
De Zeeuw 1985), (q;,/3, 7) = (2,4,0), which is the sole (known) ellipsoidal 
density profile where a test particle admits three global integrals of motion. 
Finally, it includes the Hernquist (1990) density profile, (a,/?, 7) = (1,4,1), 
which closely approximates the de Vaucouleurs r^/^ law for elliptical galaxies. 
In dealing with the formation of dark matter haloes from hierarchical clus- 
tering in both CDM and ACDM scenarios, recent high-resolution simulations 
allow (a,/3,7) = (3/2,3,3/2), hereafter quoted as MOA density profile, as a 
best fit (e.g., Ghigna et al. 2000; FMOl; KLAOl; FM03), as first advocated 
by Moore et al. (1998) and MOA99^. But for a different point of view, 
concerning the trend near the centre of the system, see e.g. Miicket & Hoeft 
(2003). 

In addition, purely dark matter structures which fulfill Jeans equation, 
exhibit —3 < 7 < — 1 for density profiles following an exact power-law, 
p oc r"'*', and this constraint weakens only slightly for a more general mass 
distribution where the density-power slope, 7(r), is a function of the radius; 
if otherwise, the system cannot be considered as in equilibrium and/or the 
effects of baryionic component have to be investigated (Hansen 2004). 

Leaving aside peculiar situations such as the occurrence of major mergers, 
the average evolution may be approximated as self-similar to a good extent. 
Accordingly, a single halo may be characterized by two parameters: the 
(fiducial) total mass, M, and a dimensionless quantity, 6, related to the 
amplitude of the density perturbation at the collapse (NFW97, FMOl). The 

-•^More precisely, a slope a — 1.4 was derived by Moore et al. (1998), while the value 
a — 1.5 was established in MOA99. 



scaling density, po, and the scaling radius, ro, may also be expressed in terms 
of M and 6 (NFW97, FMOl). 

Though Eq. (1) implies null density at infinite radius, profile fits are nec- 
essarily performed within the virialized region of a halo, bounded by a trun- 
cation radius. The presence of neighbouring systems makes the tidal radius 
as an upper limit. On the other hand, isolated objects cannot extend outside 
the Hubble sphere of equal mass. The region enclosed within the trunca- 
tion boundary has to be intended as representative of the quasi static halo 
interior, leaving aside the surrounding material which is still infalling. Nu- 
merical simulations show that the quasi static halo interior is defined by a 
mean density, paoo ~ 200pcrit (e.g.. Cole & Lacey 1996; NFW97; FMOl), 
where pcrit is the critical density of the universe. An alternative definition is 
found in KLAOl, where the quasi static halo interior has same mean density 
as predicted by the top-hat model. 

Given a set of simulated, dark matter haloes, the choice of a fitting density 
profile, expressed by Eq. (1), implies the following steps (e.g., Dubinski & 
Carlberg 1991; FMOl; KLAOl; FM03). 

(a) Select a choice of exponents (q;,/9,7), for defining the universal density 

profile. 

(b) Use a nonlinear least-squares method to determine the best fit for the 

scaling density, po, and the scaling radius, tq, with regard to each 
simulation. 

(c) Determine the scaled, averaged density profile, and related values of the 

scaling parameters, (Tq, Pq). 

(d) Particularize the fitting formula, expressed by Eq. (1), to (ro,po) = 

(ro, Po), and calculate the parameters of interest, including the residuals 
related to the scaled, averaged and fitting density profile. 

At present, no general consensus exists on the details of the above men- 
tioned procedure. The scaling density and the scaling radius may be con- 
strained to yield M{R) = Mtm, where M{R) is the mass within the trun- 
cation radius, R, related to the fitting density profile, and Mtm is the mass 
within the virialized region (e.g., FM03). An average all over the simulations 
may be performed with regard to the scaling density and the scaling radius, or 



any other two equivalent parameters (e.g., FMOl). The minimization proce- 
dure may be apphed on the sum of either the squares of logarithmic residuals, 
[log(psim/Po) -log(p«m/po)]^ (e.g., Dubiuski & Carlberg 1991; FM03), or the 
absolute values of logarithmic residuals, | log(psjm/po) ~" log(PMm/Po)| (e-g-; 
KLAOl), where psim and puni denote simulated and universal density profile, 
respectively, related to an assigned radius, r. For further details on the fitting 
procedures see e.g., Fukushige et al. (2004), Tasitsiomi et al. (2004). 

Though the above mentioned method allows a selection between different, 
universal density profiles (e.g., FMOl, FM03), it does not seem to hold in 
general (e.g., KLAOl). In other words a simulated density profile may be 
fitted, to an acceptable extent, by universal density profiles with several 
choices of exponents, (a, /9, 7), appearing in Eq. (1). Further investigation on 
additional criterions, in fitting universal to simulated density profiles, could 
be useful in lowering the above mentioned degeneracy. 

Data from a statistical sample of about five thousands of simulated, dark 
matter haloes (Bullock et al. 2001), are consistent with a lognormal dis- 
tribution of the concentration, "^tm = R/f^o, i-e. the ratio of the truncation 
radius, R, to the scaling radius, ro, with regard to NFW density profiles. The 
distribution is related to masses within a range (0.5-1.0) xlO"/i~^M0, where 
11 < n < 14 and n is an integer. The scatter is large, of about fXiogs^^^ = 0.18 
(Bullock et al. 2001). 

The existence of a lognormal distribution is a necessary, but not sufficient 
condition, for the validity of the central limit theorem. In this view, the 
concentration is related to the final properties of a simulated halo, which are 
connected with the initial conditions, ai, a2, ■■■, ««> by a transformation, 
Sirra = Q;ia2---ttn, as in dealing with the process of star formation, where the 
stellar mass follows a lognormal distribution (for further details, see Adams 
& Fatuzzo 1996; Padoan et al. 1997). 

If, for a selected mass range, the density profile is assumed to be universal, 
then the scaled physical parameters, related to e.g., mass, moment of inertia, 
and potential energy, depend on the scaled radius only (e.g., Caimmi & 
Marmo 2003). Accordingly, the distribution of a scaled physical parameter, 
to a first extent, is expected to be normal or lognormal. 

There is a well known analogon of the above procedure in the theory of 
errors. Let a physical quantity, S, be directly measured, and then follow 
a normal distribution, characterized by an expected value, S*, and a rms 
error, 0"=. Let one other physical quantity, \E', depend on the former one, 
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\1/ = 0(H), and then be indirectly measured. As a result of the theory of 
errors, the physical quantity, \E', at least to a first extent, also follows a 
normal distribution, characterized by an expected value, \E'* = 0(2*), and a 
rms error, aq, = |((90/(9S)=. |o"h. 

The current investigation is aimed to provide an additional criterion in 
fitting universal to simulated density profiles. The related procedure starts 
from a set of high-resolution simulations, where (i) both mass and radius 
of the virialized region are known for each sample halo; (2) scaled density 
profiles are averaged over the whole sample; (3) scaling parameters, (ro,po), 
are deduced with regard to different choices of exponents, {a,(3,'~f), which 
appear in the fitting formula, expressed by Eq. (1). 

A fitting, scaled density profile, is defined by a choice of exponents, 
{a,(3,'~f), and a scaled truncation radius, S, which allows the calculation 
of the remaining scaled parameters. Under the assumption of invariant, 
scaled truncation radius, a fitting, scaled density profile represents a fam- 
ily of infinite density profiles, hereafter quoted as fitting scaled halo, each 
related to a particular choice of scaling parameters, (ro,po)- Accordingly, 
the generic member of the above mentioned family is defined by three expo- 
nents, (q;,/5, 7), a scaled truncation radius, H, and two scaling parameters, 

(^o,Po)- 

A similar situation occurs for polytropes (e.g., Caimmi 1980), where a 
scaled density profile represents a family of infinite density profiles. It de- 
pends on one exponent, n (the polytropic index), and a scaled radius, S 
(where the density falls to zero). The generic member of the family depends, 
in addition, on two scaling parameters, (-R, A), which represent the radius 
and the central density, respectively. 

Given a sample of simulated, dark matter haloes, let us define a mean, 
scaled density profile, as the result of some averaging procedure over the 
whole sample. Then a mean, scaled density profile represents its parent set 
of simulations and, in the following, it shall be quoted as mean scaled halo. 

With regard to a selected, fitting, scaled density profile, (q;,/3,7), the 
dependence of the scaling parameters, (ro,po), on a pair of independent pa- 
rameters, (M, (5), can be deduced from the mean, scaled halo (e.g., FMOl). 
Then a scaled mass, M/Mq, can be explicitly expressed and compared with 
its counterpart related to the fitting, scaled density profile. It, in turn, allows 
the calculation of the scaled truncation radius, S, and other scaled parame- 
ters, which define the fitting scaled halo. 
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On the other hand, in deahng with a generic, simulated halo, both mass 
and radius of the virialized region, Mtm and Rtm-, are known as computer 
outputs. Accordingly, a scaled radius, ^tm = f'tmlf'^-, and a scaled mass, 
Mtrn/^Qi can be determined together with additional scaled parameters, and 
the deviations from their counterparts, related to the fitting, scaled halo, may 
be analysed. In particular, the arithmetic mean, rj, the standard deviation 
from the mean, agrj-, and the standard deviation from the standard deviation 
from the mean, a^jj, may be calculated for a selected scaled parameter, 77. 

The fitting density profile to a fixed simulated halo, is defined by three 
exponents, (q;,/5,7), a scaled truncation radius, S,trn, and two scaling pa- 
rameters, (ro,po)- It shall be hereafter quoted as fitting halo. Unlike the 
fitting, scaled halo, fitting haloes exhibit different scaled truncation radii, 
related to different simulated haloes, while the remaining parameters are left 
unchanged. 

At this stage, it is possible to see to what extent different, fitting, scaled 
density profiles, (a,/?, 7), make scaled parameters, r^, related to each simu- 
lated halo, deviate from their counterparts, r/*, related to the fitting, scaled 
halo. In other words, one is able to recognize if the inequality, fj — asri < 1]* < 
fj+crsrj, is fulfilled. Finally, the best fitting density profile among the ones un- 
der consideration, is chosen as minimizing the ratio, \xjj\/asrf = \v ~''l*\/^sri- 

The current investigation shall be limited to samples of recent, high- 
resolution, virialized structures where (a) the sample is homogeneous i.e. 
related to a fixed cosmological model; (b) the sample is not extremely poor 
i.e. the number of objects exceeds ten; (c) the values of the scaling density, 
scaling radius, virial mass, and virial radius, are reported or may be deduced 
from the results. 

The above conditions are satisfied by two samples, each made of a dozen 
of runs, namely FMOl and KLAOl. With regard to the latter, the twelve 
runs studied therein are in fact three sets of simulations of only four dark 
matter haloes with resolution varied in each set. They cannot be treated as 
twelve independent runs but, on the other hand, they can be conceived as 
measures of a same physical quantity, but using different methods. 

The two sets of runs differ in many respects, namely: 1) cosmological 
model; 2) criterium in making subsets of runs; 3) criterium in ending simu- 
lations; 4) mass range; 5) definition of the virial radius; 6) scaling between 
NEW and MOA density profiles; 7) choice of the pair of independent param- 
eters, i.e. (M, (5) or {po,ro). For further details, see FMOl and KLAOl. 
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Both NFW and MOA density profiles are fitted to simulated density 
profiles for the samples under consideration (FMOl, KLAOl). Accordingly, 
the above mentioned criterion shall be used in the present paper, to see what 
is the best fit, among NFW and MOA, to each set of simulations. 

The main limit of the current approach lies in the assumption of scaled 
density profiles, related to an invariant, scaled (truncation) radius, S, and 
other scaled parameters depending only on S. In general, the scaled radius, 
which has the same formal definition as the concentration (e.g., NFW97), 
depends on both the mass and the redshift. More precisely, the concen- 
tration is lowered for increasing mass (constant redshift) and redshift (con- 
stant mass), with a milder/steeper dependence for CDM/ACDM cosmologi- 
cal models (Bullok et al. 2001). An investigation based on a large statistical 
sample (Bullok et al. 2001) has shown that, within a ACDM scenario, the 
intrinsic spread in concentration related to a mass bin of distinct haloes, is 
comparable to the systematic change in the mean value of concentration re- 
lated to the above mentioned mass bin, across the entire mass range studied 
therein (10^^ < M/Mq < 10^^). 

It will be shown that, for both FMOl and KLAOl simulations, the intrinsic 
spread in concentration is dominant over the systematic change in the mean 
value of concentration on a mass bin, across the entire mass range studied 
therein. Accordingly, the fitting density profile may be considered, to an 
acceptable extent, as related to an invariant scaled radius. In general, it is 
the case for a sufficiently narrow mass range. 

The current paper is organized in the following way. Useful formulae 
related to NFW and MOA density profiles are summarized in section 2. The 
fitting, scaled haloes, related to FMOl and KLAOl set of simulations, with 
regard to both NFW and MOA density profiles, are determined in section 
3. The deviations of simulated haloes from their fitting counterparts, in 
connection with a number of scaled parameters, is also analysed therein. The 
following section 4 is dedicated to a discussion, within which some features of 
the early evolution of fitting haloes are discussed, in the limit of the spherical 
top-hat model. Some concluding remarks are drawn in section 5. Further 
investigation on a few special arguments is performed in the Appendix. 
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2 NFW and MOA density profiles 

With regard to the family of density profiles, expressed by Eq. (1), let us 
define a scaled density, /, and a scaled radius, C,, as: 

/(i) = i ; (2) 

(3) 

where the normalization, /(I) = 1, makes po and tq be the density and the 
radius (i.e. radial coordinate), respectively, of a reference isopycnic surface, 
and S corresponds to the truncation isopycnic surface, or the truncation 
radius, R. On the other hand, the normalization currently used in the lit- 
erature takes Pq = 2^po, in particular {p'q)nfw = 4po and {p'qJmoa = 2po- 
The choice of exponents, (a,/?, 7) = (1,3,1), (3/2,3,3/2), selects NFW and 
MOA density profiles, respectively, from Eq. (1). 

The explicit expression of a number of scaled parameters, related to global 
or local properties of the parent density profile, are listed in Tab. 1. Local 
properties depend on the scaled radius, C,, and global properties depend on 
the scaled truncation radius, S, which has the same formal definition as the 
concentration (e.g., NFW97). 

3 Mean and fitting dark matter haloes 

It can be shown that both NFW and MOA density profiles provide an ac- 
ceptable fit to high-resolution simulations (e.g., FMOl; KLAOl), with the 
possible exception of scales of the order of cluster of galaxies, where MOA 
density profiles seem to be preferred with respect to NFW density profiles 
(e.g., FM03). 

The explicit expression of the scaling density, po, as a function of the 
independent parameters, {M,6), prescribed in FMOl, by averaging over the 
whole set of simulations, reads: 



7k^ 7 

{Cp)nfw,fm = -j^ ; {Cp)moa,fm = — ; (4b) 
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Table 1: Comparison between functions (local properties) and profile pa- 
rameters (global properties), related to NFW and MOA density profiles, 
respectively. The profile parameters depend on a single unknown, i.e. the 
scaled radius, S. The profile parameter, z/j, is related to the special case of 
constant rotational velocity on the equatorial plane. Rigidly rotating con- 
figurations correspond to uj = uj. Caption of symbols: M - total mass 
within the truncation isopycnic surface; Mq - mass of a homogeneous region 
with same density and boundary as the reference isopycnic surface; p - mean 
density within the truncation isopycnic surface; Veq(S), (fo)eg - rotational 
velocity with respect to the centre of mass, at a point placed on the trunca- 
tion and reference isopycnic surface, respectively; / - moment of inertia; R - 
radius; Egei - self potential-energy; G - constant of gravitation; J - angular 
momentum. 
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and the analogon for the scahng radius, tq, reads: 

/ /\# \ 2/3 

r, = CrS-'/'(^-^j kpc ; (5a) 

{Cr)NFW,FM = 2 10" '" k^ ] {Cr)MOA,FM = 2 10" ; (5b) 

where the normahzation constant, ki, provides a connection between NFW 
and MOA density profiles, in fitting the results of simulations. For further 
details, see FMOl and Caimmi & Marmo (2003, Appendix B). 

With regard to a selected, scaled density profile (NFW or MOA), scaled 
parameters related to simulated haloes may be calculated via the scaling 
parameters, {rQ,pQ), defined by Eqs. (4) and (5) which, in turn, correspond 
to fixed choices of independent parameters, {M,6). For sake of brevity, let 
us define any parameter, related to a fitting, halo, as fitting parameter (e.g., 
the mass of a fitting halo is referred to as fitting mass). 

The value of the dimensionless parameter, 6, is considered to reflect an 
amplitude of the density perturbation at turnaround and, for this reason, 
it can be thought of as constant during the evolution of a halo (e.g.. Cole 
& Lacey 1996; NFW97; FMOl). From the standpoint of top-hat, spherical 
density perturbation, it is related to both the mass and the peak height, as 
shown in Appendix A. 

The combination of Eqs. (4) and (5) yields: 

MSI 



r 



where Mq is the mass of a homogeneous region, with same density and bound- 
ary as the reference isopycnic surface, (ro, po)- The shape factor, z/m, depends 
on the scaled radius, S, as shown in Tab. 1. Then the last may be determined, 
with regard to a selected, fitting density profile. 
Let us define a dimensionless parameter, /t, as: 

where R = Sro is the radius of the truncation isopycnic surface. The combi- 
nation of Eqs. (5) and (7) yields: 
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which, in turn, is related to the whole set of simulations, and depends on the 
fitting scaled radius, S. 

The simulated parameters e.g., mass, mean density, radius, of the viri- 
alized configuration, and scaled radius, shall be labeled by the index, trn, 
where trn = 200, vir, according if FMOl or KLAOl runs are involved. This is 
why different definitions of radius of the virialized configuration, have been 
used in FMOl and KLAOl. 

The scaling density, po, and the scaling radius, tq, are taken as fundamen- 
tal in KLAOl, and then no counterpart to Eqs. (4) and (5) is provided therein. 
On the other hand, the results of numerical simulations (e.g., NFW97; FMOl; 
KLAOl; Bullock et al. 2001; FM03) provide additional support to the idea, 
that density profiles of dark matter haloes in hierarchically clustering uni- 
verses have the same shape, independent of the halo mass, the initial density 
perturbation spectrum, and the values of cosmological parameters. Accord- 
ingly, we suppose that Eqs. (4) and (5) hold even in averaging scaled density 
profiles from KLAOl simulations, but different values must be assigned to 
the coefficients, Cp and Cr- In doing this, the procedure will depend on the 
density profile (NFW or MOA) under consideration. 

With regard to NFW density profiles, the scaled radius, C,mr = rmr/ro, 
is provided for each run in KLAOl, and the mean density of the virialized 
configuration, p^r = 3M^,jr/(47rr^j^), together with the scaling radius, vq, may 
be deduced from KLAOl results. The twelve runs from KLAOl correspond 
to virial masses of the same order, which allows a comparison with four runs 
executed by FMOl i.e. 4M0, 2M0, 2M1, 2M2, where the virial masses are 
also of the same order. 

It can be seen that the related, averaged, scaling radius is: 

{ro)NFW,KLA = 25.675 kpc ; {ro)NFW,FM = 12.87 kpc ; (9) 

and the ratio equals two within the uncertainty of the results. Then we 
assume that the value of the constant, Cr, appearing in Eqs. (5), doubles its 
counterpart related to FMOl simulations, that is: 

iCr)NFW,KLA = ^lO-'^^K' ; (Cr) MOA,KLA = ^ IQ-^^' • (10) 

The parameters, M^r, f^vir, and 6, are intrinsic to simulations (e.g., FMOl) 
and for this reason do not depend on the fitting density profile. Accordingly, 
the combination of Eqs. (5a) and (10) yields: 

{ro)MOA = ki{ro)NFW ] (H) 
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where ki = 2.275 according to Caimmi & Marmo (2003, Appendix B). 

The fitting scaled radius, S, is determined by averaging the results from 
KLAOl runs, as: 

S = (e.^r)jvFH^ = 13.50833 ; (12) 

in the mass range (0.68 — 2.10) x lO^'^h^^M.Q, which is consistent with the 
expected value of the lognormal distribution deduced from a statistical sam- 
ple of about two thousands of dark matter haloes in the mass range (0.5 — 
1.0) X 10^'^h-^MQ (Bullock et al. 2001). 

Then the combination of Eqs. (6) and (10) produces: 

iCp)NFW,KLA = -r- {^M)NFw{Cr)NFW,KLA 5 (1^) 

where M = umMq according to the results listed in Table 1. The combination 
of Eqs. (12) and (13) yields {Cp)nfw,kla = 0.209911. Replacing the trun- 
cation radius with the virial radius, the following relation is deduced from 
Table 1 : 

ln(l + ^^ir' 



M„v = 12Mn 



which allows the calculation of the scaling mass, Mq, the fitting mass, M, 
and then the remaining parameters, for each run in connection with NFW 
density profiles. 

In dealing with MOA density profiles, the constant, Cr, and the scaling 
radius, tq, are expressed by Eqs. (10) and (11), respectively, and the constant, 
Cp, takes the expression: 

2 75 

{Cp)mOA,KLA = -j-^{Cp)nFW,KLA = 7T^(^M)Ari?VK 5 i^^) 

where the profile parameter, um, rnay be calculated using the results listed 
in Table 1 together with Eq. (12), yielding {Cp)moa,kla = 0.0356551. The 
combination of Eqs. (6), (10), and (15) produces: 

{^m)mOA = -il^M)NFW ; (16) 

which implies {Mq)moa = '^{Mq)j^piv provided the fitting mass, M, is kept 
fixed passing from MOA to NFW density profiles and vice versa. To this 
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respect, it has already been said that the parameters, M^^, Vmr, and S, 
are also left unchanged. Then the remaining parameters may be calculated 
following the same procedure used in connection with FMOl simulations, as 
outlined in Appendix B. 

3.1 Fitting, scaled dark matter haloes related to MO A 
and NFW density profiles 

With regard to MOA density profiles, the combination of Eqs. (4b), (5b) and 
(6) yields: 

-M = Y4^ , (17) 

and the comparison with the explicit expression of the profile factor, um, 
listed in Table 1, produces: 



/375 



2/3 



(18) 



With regard to NFW density profiles, similarly, we have: 



375 

i^M = ^ ; (19) 

and the comparison with the explicit expression of the profile factor, um, 
listed in Table 1, yields: 

^-,„^ = l + Hi. (20) 

1 + S 1 + S 287r ^ ^ 

The knowledge of the scaled radius, S, via Eqs. (18) and (20), allows the 
calculation of the profile parameters, Up, z/gg, z//, Ugei, ^j, according to the 
explicit expressions listed in Table 1. The related physical parameters, p, 
Veq, I, Eseh J, may, in turn, be calculated, according to the definitions listed 
in Table 1 and via Eqs. (4), (5), provided the independent parameters, (M, 6), 
are assigned. 

The scaled radius, C,max, where either NFW or MOA velocity profile, 
Veq{^), related to centrifugal support along a selected radial direction, attains 
its maximum value, may be calculated by applying the standard methods of 
analysis and then solving a transcendental equation. 
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parameter 


FMOl 


KLAOl 




NFW 


MOA 


NFW 


MOA 


lHj 


9.20678 


3.80693 


13.50833 


5.43586 


<,max 


2.16258 


1.24968 


2.16258 


1.24968 


^P 


0.00417037 


0.0319486 


0.00140677 


0.0115409 


^P 


0.0218504 


0.154536 


0.00848858 


0.0651335 


T^M 


17.05231 


8.52616 


20.92379 


10.46189 


T^eq 


0.893929 


0.898766 


0.647442 


0.833159 


^I 


0.0554130 


0.0892287 


0.0498080 


0.0799089 


^sel 


1.12890 


0.561202 


1.10549 


0.632324 


T^J 


0.139180 


0.146741 


0.128641 


0.135717 


^rot 


0.333333 


0.333333 


0.333333 


0.333333 


K 


2.30269 


2.10091 


11.57498 


10.13895 



Table 2: Values of the scaled radius, S, the scaled radius where maximum 
centrifugal support along a selected, radial direction occurs, ^maxi the profile 

z//, Useh ^Ji ^rot, and the dimensionless parame- 



parameters, Up, Up, um, 



'eq) 



ter, K, related to fitting NFW and MOA density profiles, to both FMOl and 
KLAOl simulations. The profile parameters, vj and Vrot) are related to the 
special case of constant rotational velocity on the equatorial plane. Rigidly 
rotating configurations correspond to vj = uj, u^ot = ^i- The profile pa- 
rameter, Ueq, attains values which are very close each to the other, for FMOl 
simulations. 

Numerical values of the scaled radius, S, the scaled radius where maxi- 
mum centrifugal support along a selected, radial direction occurs, ^max, and 
the profile parameters, Up, Up, um, i^eq, ^i, ^pot, ^j, Kot, are listed in Table 2 
for fitting NFW and MOA density profiles, to both FMOl and KLAOl simu- 
lations. The profile parameter, u^q, appears to attain values which are very 
close each to the other, for FMOl simulations. 

3.2 Deviation of simulated dark matter haloes from 
their fitting counterparts 

The above results allow the comparison between simulated and fitting, dark 
matter haloes. Values of some relevant parameters related to simulations 
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run 


m/Mio 


^start 


^end 


^200 


M200/M10 


r2oo/kpc 


6-^ 


16M0 


3.0 10-2 


18.8 


0.0 


873170 


2.6 104 


1.7 103 


1.0 


16M1 


6.0 10-2 


18.5 


0.0 


1279383 


7.8 104 


2.4 103 


0.4 


16M2 


6.1 10-2 


20.4 


0.0 


1322351 


8.0 104 


2.4 103 


0.6 


8M0 


3.7 10-3 


22.3 


0.58 


745735 


2.8 103 


4.8 102 


2.5 


8M1 


7.6 10-3 


22.2 


0.63 


1186162 


9.0 103 


7.2 102 


1.0 


8M2 


7.6 10-3 


23.9 


0.59 


1015454 


7.7 103 


7.0 102 


3.0 


4M0 


4.7 10"^ 


25.9 


1.6 


559563 


2.7 102 


1.3 102 


10.0 


4M1 


9.5 10-4 


25.9 


1.6 


846301 


8.0 102 


2.0 102 


3.0 


4M2 


9.5 10-4 


27.4 


1.2 


697504 


6.6 102 


2.2 102 


6.0 


2M0 


5.9 10-5 


29.7 


2.1 


643151 


6.6 10^ 


6.2 10^ 


35.0 


2M1 


1.2 10-4 


29.7 


2.2 


957365 


1.1 102 


8.5 10^ 


12.0 


2M2 


1.2 10-4 


30.9 


1.8 


923545 


1.0 102 


9.6 10^ 


30.0 



Table 3: Values of some parameters related to simulated, dark matter haloes, 
according to FMOl, for a standard CDM model with Hq = 50 km s~^ Mpc-^, 
Q = 1, and (Xg = 0.7. The particle masses are equal, and the total number 
of particles for each simulation is (2.0 — 2.1) x 10^. Captions: m - mass of a 
single particle; z - redshift (at the start and end of simulation); iV2oo - total 
number of particles within the sphere where p = 200pcrit {Pcnt is the critical 
density); M200 - mass enclosed within the above mentioned sphere; r2oo - 
radius of the above mentioned sphere; 5 - dimensionless parameter related to 
the amplitude of the density perturbation at the collapse. The mass unit is 
Min = 10i°M„. 



with high resolution, performed by FMOl, are listed in Table 3. 

The framework is a standard CDM model with Hq = 50 km s^^ Mpc-^, 
Q = 1, and as = 0.7. The particle masses are equal, and the total number 
of particles for each simulation is (2.0 — 2.1) x 10^. For further details, see 
FMOl. The ending redshift, z^nd > 0, is determined so that the truncation 
outside the sphere did not influence the proflle around r2oo- Then the data 
listed in Table 3 do not make a homogeneous set, as simulations related to 
Zend = satisfy a different condition in respect of Zend > 0. To this aim, 
computations related to Zend = would have been continued until the above 
mentioned condition would be satisfled at some Zend < in the future. 

Values of some relevant parameters related to simulations with high res- 
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run 


m 
Mio 


Zend 


N ■ 

1 y vtr 


V Mio ) 


\ Mio J 


M„ir 

Mio 


rvir 

kpc 


? 


Al 


2.3 10"^ 





1.2 10^ 


294 


259 


286 


367 


Y 


A2 


1.9 10-2 





1.5 10^ 


302 


268 


300 


373 


Y 


A3 


1.6 10-1 





1.9 10^ 


322 


287 


286 


366 


N 


Bl 


1.7 10-^ 





1.0 10^ 


184 


157 


171 


307 


Y 


B2 


1.7 10-^ 





1.5 10^ 


2.71 


2.39 


157 


304 


N 


B3 


1.7 10-^ 


1 


7.1 10^ 


125 


116 


121 


344 


Y 


CI 


1.7 10-^ 





1.1 10^ 


201 


173 


186 


321 


Y 


C2 


1.1 10-2 





1.6 10^ 


190 


163 


171 


314 


Y 


C3 


1.7 10-^ 


1 


5.0 10^ 


88 


82 


97 


297 


N 


Dl 


1.7 10-^ 





1.3 10^ 


236 


206 


214 


336 


Y 


D2 


1.1 10-2 





2.0 10^ 


236 


205 


214 


334 


Y 


D3 


1.7 10-^ 


1 


7.9 10^ 


139 


129 


137 


350 


Y 



Table 4: Values of some parameters related to simulated, dark matter 
haloes, according to KLAOl, for a ACDM model with Hq = 70 km s~^ 
Mpc~i, Oq = 1 — Aq = 0.3, and ag = 0.9. All simulations were started at 
Zstart = 60. Captions: m - mass of a single particle; Zend - redshift at the 
end of simulation; N^r - total number of particles within the sphere where 
P = Pcrit^o^TH {pcrit is the Critical density and 6th is the density excess 
predicted by the top-hat model); M^ir - total mass within the above men- 
tioned sphere; r^i^ - radius of the above mentioned sphere; {Nmr^TT-f^io)^ - 
upper and lower value of N^irfn/MiQ deduced from the data. The mass unit 
is Mio = 10"'^° Mq. a positive answer to the question mark means that the 
inequality, Mj^^ < M^r < M^^^ is satisfied. The related, explicit expression, 
is shown by Eqs. (22). 
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olution, performed by KLAOl, are listed in Table 4. The framework is a 
ACDM model with Hq = 70 km s"^ Mpc"\ I^q = I-Aq = 0.3, and ag = 0.9. 
All simulations were started at Zgtart = 60. The twelve runs are, in fact, 
three sets of simulations of only four haloes with resolution varied in each 
set. Thought they cannot be considered as twelve independent runs, still they 
may be conceived as measures of a same physical quantity, using different 
methods. For further details, see KLAOl. 

It is worth of note that the virial radius in KLAOl is not defined as the 
radius, r2ooj within which the mean density is 200 times the critical density, 
as done by e.g., NFW97 and FMOl. On the other hand, the virial radius is 
intended therein as the radius, Vmr, within which the mean density is equal 
to the density predicted by the top-hat model, pmr = STH^oPcrit, where Sth 
is the density excess predicted by the top-hat model and pcrit is the critical 
density. In the case of Qq = 0.3 cosmologies, it can be seen that r^ir ~ 1.3r2oo 
(KLAOl). 

In general, the scaled radius, H (or concentration with regard to NFW 
density profiles), is lowered for increasing mass (constant redshift) and red- 
shift (constant mass). The mass range is large (6 ■ 10^^ < M/Mq < 8 ■ lO^'') 
for FMOl simulations, but the use of a CDM cosmological model makes the 
concentration mildly depend on the mass (at constant redshift). On the other 
hand, the concentration is decreased for low-mass haloes, which virialize ear- 
lier, and increased for high-mass haloes, which virialize later. The net effect 
is an even milder dependence of the concentration on the mass. 

In any case, it can be said that the intrinsic spread in concentration is 
dominant over the systematic change in the mean value of the concentration 
on a mass bin (e.g., Bullock et al. 2001), across the entire mass ranges 
covered by FMOl and KLAOl simulations. Accordingly, the fitting, scaled 
density profile may be considered, to an acceptable extent, as related to an 
invariant scaled radius, S. 

The following inequalities may be useful for testing the intrinsic spread 
of some data listed in Table 3: 

^200 < M200 < M+Q ; P200 < 20(}pcrit{zend) < ptoo ! (21a) 

P200 ^ P200 < P200 5 ^^200 ^ '^soo < '^200 5 (21b) 

^200 = (^200 T ^N2oo) {m T Am) ; (21c) 



P200 — P200 



1 T AM200/M 



200 



^1 ± Ar2oo 7^200)^ 



(21d) 



,^,,r-^(-|)- ., „,, 



^200 

3 M200 

3 

200 

Pe„t(^) = 0.691785 10-^(1 + 2)3Miokpc-=^ ; (21g) 

3/2 , ,^ X -1 



P200 = -. — 3— ; (21f) 

Air To 



rl/2 M200 \ /^200\ /01U^ 



where z is the redshift, Mio = lO^^M© and, in general, Arj = 5 ■ 10""""'^ is the 
uncertainty assumed for rj = u 10~", < -u < 10. Upper and lower values are 
listed in Table 5. A positive answer to a question mark therein, with regard 
to the parameter on the left of the column under consideration, means that 
the related inequalities, among (21a) and (21b), are satisfied, and vice versa. 
Similarly, for testing the intrinsic spread of masses listed in Table 4: 

M-,<M,ir<M+^ ; (22a) 

M^,, = {N,,r T AN,ir){m T Am) ; (22b) 

where, in general, Arj = 5- 10^"^^ is the uncertainty assumed for r^ = -u 10^", 
< M < 10. The upper and lower values are also listed in Table 4. A 
positive answer to the question mark therein means that the inequality (22a) 
is satisfied, and vice versa. 

All the data produce an acceptable intrinsic spread, with the exception 
of run 2M0, in FMOl, and run B2, in KLAOl, which exhibit a substantial 
inconsistency in any case. The larger discrepancies can be due to nothing 
but printing errors. Accordingly, the value M200/M10 = 38, also deduced 
from FMOl, is assigned instead of 66, to run 2M0, and the value m/Mio = 
1.1 10~^, also inferred from KLAOl, is assigned instead of 1.7 10~^, to run 
B2. 

With regard to NEW density profiles, the scaled radius, C,mr = rmr/ro, 
is provided for each run in KLAOl, and the mean density of the virialized 
configuration, p^r = 3M„jj./(47rr^j^), together with the scaling radius, vq, may 
be deduced from the results of Table 4. The above mentioned parameters 
are listed in Table 6, together with two dimensionless parameters, 6 and K^r, 
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run 


M200 

Mio 


( N2mm\'^ 
V Mio ) 


( N2QQm\ 
V Mio ) 


? 


/t200 




'^200 


'^200 


? 


16M0 


26000 


26632 


25758 


Y 


2.70 




2.92 


2.46 


Y 


16M1 


78000 


77403 


76123 


N 


2.38 




2.69 


2.19 


Y 


16M2 


80000 


81325 


80002 


N 


1.90 




2.04 


1.74 


Y 


8M0 


2800 


2796 


2722 


N 


2.37 




2.51 


2.32 


Y 


8M1 


9000 


9074 


8955 


Y 


2.15 




2.24 


2.06 


Y 


8M2 


7700 


7768 


7667 


Y 


1.39 




1.42 


1.37 


Y 


4M0 


270 


266 


260 


N 


1.74 




1.91 


1.66 


Y 


4M1 


800 


808 


800 


Y 


2.04 




2.14 


1.93 


Y 


4M2 


660 


666 


659 


Y 


2.02 




2.10 


1.92 


Y 


2M0 


66 


38 


38 


N 


1.25 




2.20 


2.14 


N 


2M1 


110 


120 


110 


Y 


2.06 




2.08 


1.86 


Y 


2M2 


100 


115 


106 


N 


1.72 




1.63 


1.48 


N 




200pcrit{Zend) 

Miokpc"^ 


P200 


P'nm 


? 


P200 


T 


? 








Miokpc^^ 


Miokpc"-^ 


Miokpc^' 




16M0 


1.38 10-6 


1.41 10-6 


1.15 10-6 


Y 


1.26 10- 


-6 


Y 






16M1 


1.38 10-6 


1.42 10-6 


1.24 10-6 


Y 


1.35 10- 


-6 


Y 






16M2 


1.38 10-6 


1.50 10-6 


1.30 10-6 


Y 


1.38 10- 


-6 


Y 






8M0 


5.46 10-6 


6.23 10-6 


5.70 10-6 


N 


6.04 10- 


-6 


Y 






8M1 


5.99 10-6 


5.93 10-6 


5.61 10-6 


N 


5.76 10- 


-6 


Y 






8M2 


5.56 10-6 


5.52 10-6 


5.22 10-6 


N 


5.36 10- 


-6 


Y 






4M0 


2.43 10-5 


3.25 10-5 


2.52 10-5 


N 


2.93 10- 


-5 


Y 






AMI 


2.43 10-5 


2.60 10-5 


2.22 10-5 


Y 


2.39 10- 


-5 


Y 






4M2 


1.47 10-5 


1.60 10-5 


1.38 10-5 


Y 


1.48 10" 


-5 


Y 






2M0 


4.12 10-5 


3.90 10-5 


3.72 10-5 


N 


6.61 10- 


-5 


N 






2M1 


4.53 10-5 


4.75 10-5 


4.20 10-5 


Y 


4.28 10- 


-5 


Y 






2M2 


3.04 10-5 


3.15 10-5 


2.82 10-5 


Y 


2.70 10- 


-5 


N 







Table 5: Comparison between (i) the value of M200 and the upper and lower 
value of N2oom, deduced from the data listed in Table 3; (ii) the value of 
K200 deduced from the data listed in Table 3 and the related upper and 
lower values, using Eqs. (21e) and (21h), respectively; and (iii) the value 
of 200pcrit{zend) Calculated using Eq. (21g), the upper and lower value of 
P200 = 3A^2oo"^/(47rr|Qo), and the value of P200 = 3M2oo/(4vrr|oo), deduced 
from the data listed in Table 3. A positive answer to the question mark 
means that the related inequalities, among Eqs. (21a) and (21c), are satisfied 
for the values listed on the same lines of the corresponding columns. The 
related, explicit expressions, are shown by Eqs. (21c)-(21h) . 
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i^vir)NFW 


(i"o)nfw 
kpc 


Pvir 


6-' 






Miokpc"^ 


f^vtr 


Al 


17.4 


21.1 


1.38 10-6 


2.68 


16.9 


Al 


16.0 


23.3 


1.38 10-6 


3.05 


14.9 


A3 


16.6 


22.0 


1.39 10-6 


2.94 


15.8 


B\ 


15.6 


19.7 


1.41 10-6 


5.51 


14.4 


B2 


16.5 


18.4 


1.33 10-6 


5.66 


15.6 


B3 


12.3 


28.0 


7.10 10-^ 


25.0 


10.1 


CI 


11.2 


28.7 


1.33 10-6 


10.3 


8.74 


C2 


9.8 


32.0 


1.32 10-6 


14.7 


7.15 


C3 


11.9 


25.0 


8.84 10-^ 


26.7 


9.57 


Dl 


11.9 


28.2 


1.35 10-6 


7.95 


9.57 


D2 


13.4 


24.9 


1.37 10-6 


6.17 


11.4 


D3 


9.5 


36.8 


7.64 10-^ 


33.7 


6.83 



Table 6: The scaled, virialized radius related to NFW density profiles, 
{^vir)NFW, the scaling radius related to NFW density profiles, {ro)NFW, the 
mean density inside the virialized configuration, p^r, and the dimensionless 
parameters, 6~^ and K,mr, taken from twelve runs in KLAOl or deduced from 

Eqs. (5), (7), and (8). 
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which are deduced from Eqs. (5), (7), and (8). It is apparent that the twelve 
runs correspond to virial masses of the same order. 

At this stage, the deviation of simulated, dark matter haloes, from their 
fitting counterparts, may be analysed along the following steps. 

(i) Select a fitting halo among NFW and MOA density profiles. 

(ii) Select a simulated halo among the twelve runs in Tables 3 and 4. 

(iii) Calculate the fitting mass, M, using Eqs. (4), (5), (6), and the related 
values of Mtrn and rtm-, trn = 200, vir, listed in Tables 3 and 4. For a 
formal derivation, see Appendix B. 

(iv) Calculate the values of the scaling density, po, the scaling radius, tq, 
and the scaling mass, Mq. It is worth remembering the latter is the 
mass of a homogeneous region, with same density and boundary as the 
reference isopycnic surface, (ro,po)- 

(v) Calculate the scaled mass, Mtm/Mo, the scaled radius, rtm/ro, the 
scaled density, ptm/po, and the dimensionless parameter, utm, trn = 
200, vir. 

(vi) Return to (ii). 

(vii) Return to (i). 

Simulated haloes may be characterized by four scaled parameters, Mtm/Mo, 
rtrn/ro, Ptrn/po, and Ktrn- Their fitting counterparts are um, S, Up, and k,, 
respectively, which have been listed in Table 2. The fitting mass, M, the 
scaling density, radius, mass, po, rg, Mq, and the ratios of three parameters 
related to simulated haloes, to their fitting counterparts, /tj„i/'t, Mtrn/M, 
rtrn/R, are listed in Tables 7 and 8, in connection with both NFW and MOA 
density profiles, for the twelwe runs from FMOl and KLAOl, respectively. 

4 Discussion. 

The deviation of simulated dark matter haloes from their fitting counter- 
parts, can be seen in Figs. 1, 2, 3, 4, and in Figs. 5, 6, 7, 8, in connection with 
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M 
Mio 


PO 


T 


kpc 


Mio 


K2()(), 
K 


M 


'"200 




Miokpc ■ 


R 


NFW 


16M0 


2.37 10^ 


8.69 10- 


-5 


1.56 102 


1.39 103 


1.17 


1.10 


1.18 


16M1 


7.64 10^ 


6.74 10- 


-5 


2.51 102 


4.48 10^ 


1.03 


1.02 


1.04 


16M2 


9.09 10^ 


3.78 10- 


-5 


3.23 102 


5.33 10^ 


0.824 


0.880 


0.807 


8M0 


2.75 10=^ 


3.00 10- 


-4 


5.04 10^ 


1.61 102 


1.03 


1.02 


1.03 


8M1 


9.41 10=^ 


2.19 10- 


-4 


8.44 10^ 


5.52 102 


0.932 


0.957 


0.926 


8M2 


1.24 10^ 


6.11 10- 


-5 


1.37 102 


6.59 102 


0.603 


0.685 


0.554 


4M0 


3.27 102 


6.29 10- 


-4 


1.94 10^ 


1.92 10^ 


0.754 


0.825 


0.728 


4M1 


8.64 102 


7.94 10- 


-4 


2.49 10^ 


5.07 10^ 


0.886 


0.925 


0.876 


4M2 


7.19 10^ 


4.78 10- 


-4 


2.76 10^ 


4.21 10^ 


0.877 


0.918 


0.865 


2M0 


3.94 10^ 


1.49 10- 


^3 


7.18 10° 


2.31 10° 


0.943 


0.964 


0.938 


2M1 


1.18 102 


1.45 10- 


-3 


1.04 10^ 


6.94 10° 


0.893 


0.930 


0.884 


2M2 


1.22 102 


5.62 10- 


^4 


1.45 10^ 


7.17 10° 


0.746 


0.818 


0.719 


MOA 


16M0 


2.22 10^ 


1.58 10- 


-5 


3.40 102 


2.60 10^ 


1.28 


1.17 


1.31 


16M1 


7.17 10^ 


1.22 10- 


-5 


5.48 102 


8.41 10^ 


1.13 


1.09 


1.15 


16M2 


8.62 10^ 


6.76 10- 


-6 


7.09 102 


1.01 10^ 


0.903 


0.928 


0.889 


8M0 


2.58 10^ 


5.43 10- 


-5 


1.10 102 


3.03 102 


1.13 


1.09 


1.15 


8M1 


8.87 10^ 


3.55 10- 


-5 


1.85 102 


1.05 10^ 


1.02 


1.01 


1.02 


8M2 


1.10 10^ 


1.06 10- 


-5 


3.08 102 


1.30 10^ 


0.661 


0.698 


0.597 


4M0 


3.13 102 


1.12 10- 


-4 


4.28 10^ 


3.68 10^ 


0.826 


0.863 


0.798 


4M1 


8.17 102 


1.43 10- 


-4 


5.43 10^ 


9.59 10^ 


0.972 


0.980 


0.968 


4M2 


6.79 102 


8.59 10- 


-5 


6.05 10^ 


7.97 10^ 


0.961 


0.972 


0.955 


2M0 


3.71 10^ 


8.60 10- 


-6 


1.57 10^ 


4.36 10° 


1.03 


1.02 


1.04 


2M1 


1.12 102 


2.61 10- 


-4 


2.29 10^ 


1.31 10^ 


0.979 


0.985 


0.976 


2M2 


1.17 102 


9.98 10- 


-5 


3.20 10^ 


1.37 10^ 


0.817 


0.856 


0.788 



Table 7: The fitting mass, M, the scaling density, radius, mass, po, ''"o, ^o, 
and the ratios of three parameters related to simulated haloes to their fitting 
counterparts, /€20o/'^5 M200/M, r2oo/R, in connection with NFW and MOA 
density profiles, and twelve runs from FMOl. 



23 





M 
Mio 


Po 


Ma 
Mio 


K 


M 


l^vir 




Miokpc"^ 


R 


NFW 


Al 


254 


3.08 10-4 


12.1 


1.46 


1.13 


1.29 


A2 


276 


2.49 10-4 


13.2 


1.29 


1.08 


1.18 


A3 


259 


2.76 10-4 


12.4 


1.36 


1.10 


1.23 


Bl 


159 


2.39 10-4 


7.62 


1.24 


1.07 


1.15 


B2 


143 


2.60 10-4 


6.82 


1.35 


1.10 


1.22 


B3 


127 


6.62 10-5 


6.06 


0.869 


0.954 


0.910 


CI 


205 


9.93 10-5 


9.79 


0.755 


0.908 


0.829 


C2 


202 


7.02 10-5 


9.68 


0.618 


0.844 


0.725 


C3 


103 


7.59 10-5 


4.94 


0.827 


0.937 


0.881 


Dl 


228 


1.16 10-4 


10.9 


0.827 


0.937 


0.881 


D2 


215 


1.58 10-4 


10.3 


0.988 


0.996 


0.992 


D3 


165 


3.77 10-5 


7.89 


0.590 


0.830 


0.703 


MOA 


Al 


236 


5.63 10-5 


22.6 


1.48 


1.21 


1.48 


A2 


258 


4.54 10-5 


24.6 


1.36 


1.16 


1.36 


A3 


241 


5.02 10-5 


23.1 


1.41 


1.18 


1.41 


Bl 


149 


4.35 10-5 


14.2 


1.32 


1.15 


1.32 


B2 


133 


4.74 10-5 


12.7 


1.40 


1.18 


1.40 


B3 


118 


1.20 10-5 


11.3 


1.04 


1.02 


1.04 


CI 


191 


1.81 10-5 


18.3 


0.950 


0.972 


0.948 


C2 


190 


1.28 10-5 


18.1 


0.837 


0.901 


0.829 


C3 


96.6 


1.38 10-5 


9.23 


1.01 


1.00 


1.01 


Dl 


213 


2.10 10-5 


20.4 


1.01 


1.00 


1.01 


D2 


200 


2.88 10-5 


19.2 


1.13 


1.07 


1.13 


D3 


155 


6.83 10-6 


14.8 


0.812 


0.884 


0.802 



Table 8: The fitting mass, M, tlie scaling density and mass, poj and Mq, and 
the ratios of three parameters related to simulated haloes, to their fitting 
counterparts, Kmrl n, Mmr/M, Vyir/R, in connection with NFW and MOA 
density profiles, and twelve runs from KLAOl. 
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FMOl and KLAOl simulations, respectively. The following plots are repre- 
sented: scaled mass, Mtm/Mo, vs. logarithmic scaled mass, log(Mtrn/Mio); 
scaled radius, rtm/ro, vs. logarithmic scaled radius, log(rirn/kpc); scaled 
density, ptm/Po, vs. logarithmic scaled density, log[ptrn/(Miokpc~'^)]; and 
dimensionless parameter, Ktm, vs. logarithmic scaled mass, log(Mtr„/Mio); 
where tnr = 200, vir. The deviation of the above mentioned, scaled 

parameters, from their fitting counterparts, are clearly shown therein. 

The mean value, f/, the standard deviation from the mean value, agfj, and 
the standard deviation from the standard deviation from the mean value, 
asp., which are expressed as (e.g., Oliva & Terrasi 1976, Chap. V, §5.6.3): 

1 " 

v = -T.v^ ; (23) 



n 



'^ srj 



1=1 
1 1 

nn — 



1/2 



1=1 



(24) 



(^sp = —7= ; 11 = (Tsn ] (25) 

v In 



where n = 12; r] = Mtm/Mo, nm/ro, ptm/po, t^tm] tnr = 200, vir; and f], 
o'sfi, o'sfi, are shown by the vertical bars on Figs. 1-8. 

Numerical values of fj, asf^, agp, deduced by Eqs. (23), (24), and (25), 
respectively, are listed in Table 9, together with their fitting counterparts, 
rj*, which have been represented as horizontal lines in Figs. 1-8 and listed in 
Table 2, in connection with NFW and MOA density profiles. It is apparent 
that the following inequalities hold: 

f] - ua-gfj - ucsfi <1]* <f] + uasfj + uas/i ; (26) 

where m = 2, 1 for NFW and MOA density profiles, respectively, in connection 
with FMOl simulations, and the contrary holds, in connection with KLAOl 
simulations. Then the best fitting density profile (among NFW and MOA in 
the case under discussion) may be chosen as minimizing the ratio, \xjj\/(7srj = 
\r] — r]*\/(Tsrj, for the scaled parameter of interest. 

It can also be shown that a necessary condition for the detectability of 
accidental errors, Arj < a^,, is satisfied for any choice of rj listed in Table 9. 
More specifically, Ar^ and a^, represent the sensitivity error of the simulation 
and the rms error, respectively, with regard to rj. For a formal demonstration, 
see Appendix C. 
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Figure 1: Deviation of simulated dark matter haloes from their fitting coun- 
terparts, with regard to the scaled, virial mass, M200/MQ, for both NFW 
(top) and MOA (bottom) density profiles. Fitting haloes lie on the horizon- 
tal line, while simulated haloes are represented by different symbols. The 
vertical bar is centred on the mean value of plotted data, with respect to 
the ordinates (and no connection with the abscissae), and is limited by the 
standard deviation from the mean, without (inner boundary) and with (outer 
boundary) addition of about twice the standard deviation from the standard 
deviation from the mean, deduced from Eqs. (23), (24), and (25), respec- 
tively. Captions of symbols on the right correspond to FMOl runs listed in 
Tables 3, 5, and 7. 
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Figure 2: Deviation of simulated dark matter haloes from their fitting coun- 
terparts, with regard to the scaled radius, r2oo/ro, for both NFW (top) and 
MOA (bottom) self-similar, universal density profiles. Other captions as in 
Fig.l. 
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Figure 3: Deviation of simulated dark matter haloes from their fitting coun- 
terparts, with regard to the scaled density, P200/P0, ioi both NFW (top) and 
MOA (bottom) self-similar, universal density profiles. Other captions as in 
Fig.l. 
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Figure 4: Deviation of simulated dark matter haloes from their fitting coun- 
terparts, with regard to the dimensionless parameter, K200, for both NFW 
(top) and MOA (bottom) self-similar, universal density profiles. Other cap- 
tions as in Fig. 1. 
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parameter 


NFW 


MOA 


V 


V 


(^sfj 


(^sp. 


V* 


V 


O-Sjy 


O'sA 


V* 


M200/M0 


15.69 


0.538 


0.110 


17.05 


8.29 


0.309 


0.0630 


8.53 


^200/^0 


8.09 


0.443 


0.0902 


9.21 


3.69 


0.209 


0.0427 


3.81 


P200/P0 


0.0346 


0.00564 


0.00115 


0.0218 


0.195 


0.0328 


0.00669 


0.154 


'^200 


2.05 


0.100 


0.0205 


2.30 


2.05 


0.100 


0.0205 


2.10 


M„r/Mo 


20.74 


0.632 


0.129 


20.92 


11.11 


0.343 


0.0701 


10.46 


rvir/ro 


13.51 


0.807 


0.165 


13.51 


5.94 


0.355 


0.0723 


5.44 


Pvir/Po 


0.0101 


0.00155 


0.000317 


0.00849 


0.0556 


0.00855 


0.00174 


0.0651 


rvyif 


11.75 


1.03 


0.211 


11.57 


11.75 


1.03 


0.211 


10.14 



Table 9: The mean value, rj, the standard deviation from the mean value, asf^, 
and the standard deviation from the standard deviation from the mean value, 
asfi, where r] = Mtm/Mo, nm/ro, ptm/Po, i^trn, tm = 200, mr, deduced from 
simulated haloes, data from Tables 3, 5, 7 (FMOl), and 4, 6, 8 (KLAOl). The 
value of the related fitting counterparts, 1]* = um, S, Up, k, listed in Table 2, 
is also reported for comparison. The upper and lower panel are related to 
FMOl and KLAOl simulations, respectively. The left-hand and right-hand 
side are related to NFW and MOA universal density profiles, respectively. 
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It is apparent that NFW density profiles reproduce FMOl simulations to 
a lesser extent with respect to MOA density profiles, while the contrary holds 
with regard to KLAOl simulations, according to the above mentioned crite- 
rion. In any case, a better fit could be obtained by use of a different scaled 
density profile, i.e. different values of the exponents, (q;,/9,7), appearing in 
Eq.(l). 

4.1 Standard deviations deduced from the propagation 
of the errors 

Universal density profiles, expressed by Eq. (1), are currently used in fitting 
simulated density profiles (e.g., NFW; FMOl; FM03). The choice of the ex- 
ponents, (a, j3, 7), allows the scaled parameters, um = M/Mq, Up = p/po, and 
k, depend on a single scaled parameter i.e. the scaled radius, S, or concen- 
tration with regard to NFW density profiles. As the concentration exhibits a 
lognormal distribution (e.g., Bullock et al. 2001), and the above mentioned 
parameters may be considered as depending on the decimal logarithm of the 
concentration, logS, the related distribution is expected to be (at least to a 
first extent) normal, via the same procedure which leads to the propagation 
of the errors. 

The corresponding rms errors are: 

(27) 
(28) 

^logH ; (29) 

cTfc = ^kaiogs ; (30) 

where /(S) is the scaled density, expressed by Eq. (2), particularized to the 
scaled radius, ^ = S. For a formal demonstration, see Appendix D. 

A comparison between rms errors, expressed by the above relations, and 
standard deviations, listed in Table 9, needs the following steps. 

(a) Assume a rms error of the lognormal distribution of the concentration as 
in Bullock el al. (2001), aiogs = 0.18. 
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OS = ^o-iogH ; 


^^M = 3/(S)sViogH ; 


^i^p = 3 
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Figure 5: Deviation of simulated dark matter haloes from their fitting coun- 
terparts, with regard to the scaled virial mass, M„jj./Mo, for both NFW (top) 
and MOA (bottom) self-similar, universal density profiles. Fitting haloes lie 
on the horizontal line, while simulated haloes are represented by different 
symbols. The vertical bar is centred on the mean value of plotted data, with 
respect to the ordinates (and no connection with the abscissae), and is lim- 
ited by the standard deviation of the mean, without (inner boundary) and 
with (outer boundary) addition of about twice the standard deviation from 
the standard deviation from the mean, deduced from Eqs. (23), (24), and 
(25), respectively. Captions of symbols on the right correspond to KLAOl 
runs listed in Tables 4, 6, and 8. 
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Figure 6: Deviation of simulated dark matter haloes from their fitting coun- 
terparts, with regard to the scaled radius, r^,jj./ro, for both NFW (top) and 
MOA (bottom) self-similar, universal density profiles. Other captions as in 
Fig. 5. 
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Figure 7: Deviation of simulated dark matter haloes from their fitting coun- 
terparts, with regard to the scaled density, Pvir/Po, for both NFW (top) and 
MOA (bottom) self-similar, universal density profiles. Other captions as in 
Fig. 5. 
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Figure 8: Deviation of simulated dark matter haloes from their fitting coun- 
terparts, with regard to the dimensionless parameter, K^r, for both NFW 
(top) and MOA (bottom) self-similar, universal density profiles. Other cap- 
tions as in Fig. 5. 



35 



parameter 


NFW 


MOA 


V 


(Tfj 


(^fi/v* 


Xfl 


Xfj /(Tfi 


(Tfi 


(Tf^/rj* 


\X-q 


\Xfj\/(Tf^ 


M200/M0 


1.17 


0.0685 


1.36 


1.16 


0.633 


0.0742 


0.24 


0.379 


^'20o/l"o 


1.10 


0.120 


1.12 


1.02 


0.455 


0.120 


0.12 


0.264 


P200/P0 


0.0635 


0.290 


0.0128 


2.02 


0.0440 


0.285 


0.041 


0.932 


1^200 


0.413 


0.179 


0.25 


0.605 


0.377 


0.179 


0.05 


0.133 


M„r/Mo 


1.24 


0.0595 


0.18 


0.145 


0.665 


0.0636 


0.65 


0.977 


Tvir/ro 


1.62 


0.120 








0.650 


0.120 


0.5 


0.769 


Pvirl PO 


0.00254 


0.299 


0.00161 


0.634 


0.0192 


0.295 


0.0095 


0.495 


f^vir 


2.30 


0.179 


0.18 


0.0783 


1.82 


0.179 


1.61 


0.885 



Table 10: The rms errors, 0"^^, the ratios of rms error to the expected value of 
the related distribution, (Tfj/r]*, the absolute errors. 



\Xr. 



ratios of absolute error to rms error, \xfj\/afj, where r] 



- \v ^ V*\y ^^^ the 

Mtrn/Mo, r-trn/ro, 

Ptrn/Po, i^trni tm = 200, vir, deduced from simulated haloes, data from Ta- 
bles 3, 5, 7 (FMOl), and 4, 6, 8 (KLAOl). Values of parameters related to 
their fitting counterparts, f]* = um, 2, Up, k, are taken from Table 2. The 
upper and lower panel are related to FMOl and KLAOl simulations, respec- 
tively. The left-hand and right-hand side are related to NFW and MOA 
density profiles, respectively. 

(b) Assume a scaled radius related to the expected value of the lognormal 

distribution of the concentration, S* = exp^o(log'=')*5 ^^ listed in Table 
9 for each universal density profile and each set of simulations. 

(c) Divide the rms errors, a^, r] = E, um, i^p, k, by the square root of 

the number of measures that have been averaged in calculating the 
standard deviations from the mean, which is equal to vT2- 

(d) For a selected, universal density profile and a selected set of simulations, 

use the mean values, rj, and their fitting counterparts, r]*, listed in Table 
9. 

The rms errors, (Tfj, and the ratios of rms error to the expected value of the 



related distribution, afi/rj*, the absolute errors. 



Xr, 



Irj — r]*], and the ratios 



of absolute error to rms error, \xfj\/afj, are listed in Table 10 with regard to 
FMOl (upper panel) and KLAOl (lower panel) simulations, related to NFW 
(left-hand side) and MOA (right-hand side) universal density profiles. 
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The rms errors appear to be systematically larger than the related stan- 
dard deviations listed in Table 9, essentially for three orders of reasons. First, 
the statistical significance of the samples considered is low, owing to the small 
number of objects (A^ = 12). In addition, KLAOl runs make in fact three 
sets of simulations of only four dark matter haloes with resolution varied in 
each set. For this reason, they cannot be treated as twelve independent runs, 
and a lower standard deviation is expected. 

Second, the rms errors calculated by use of Eqs. (27)-(30) need a normal 
distribution for the related random variables, which holds to a good extent 
only if the fluctuations are sufficiently small, to neglect higher-order terms, 
with respect to the first, in the related series developments. For further 
details, see Appendix D. 

Third, the comparison between rms errors and standard deviations should 
be valid only in connection with the cosmological model, and the universal 
density profile, used by Bullock et al. (2001) in building up the statistical 
sample of simulated dark matter haloes, from which a lognormal distribu- 
tion of the concentration has been deduced, with rms error (T\ogEtr„ = 0.18. 
It needs Aq = 0.7, Qq = 0.3, h = 0.7, as = 1.0, and a NFW density pro- 
file. In fact, larger discrepancies occur for FMOl simulations, where Aq = 0, 
Qq = 1, h = 0.5, and as = 0.7. On the other hand, KLAOl simulations were 
performed with the same choice of cosmological parameters as in Bullock et 
al. (2001), with the exception of as = 0.9. Then we expect that a lognor- 
mal distribution of concentration for an assigned mass bin occurs for any 
plausible choice of cosmological parameters, and its expected value and rms 
error do not change dramatically for any plausible variation of cosmological 
parameters. 

The existence of a lognormal distribution is a necessary, but not sufficient 
condition, for the validity of the central limit theorem. In this view, the 
concentration is related to the final properties of a simulated halo, which are 
connected with the initial conditions, ai, 0:2, ..., ««, by a transformation, 
S = q;iq;2...q;„, as in dealing with the process of star formation, where the 
stellar mass follows a lognormal distribution (Adams & Fatuzzo 1996; Padoan 
et al. 1997). With regard to dark matter haloes within a fixed mass bin, an 
interpretation of the lognormal distribution, depending on the concentration, 
in terms of the central limit theorem, is outlined in Appendix E. 

It is worth recalling that standard deviations of scaled parameters from 
their mean values, have been deduced directly from the results of simula- 
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tions, with regard to a selected, fitting density profile. On the other hand, 
rms errors of scaled parameters have been deduced from their dependence 
on the (decimal) logarithm of the concentration, log S, and the lognormal 
distribution of the concentration. 

An inspection to Table 10 shows that the ratio of absolute error to rms 
error, is closer to zero for MOA density profiles with regard to FMOl sim- 
ulations, and for NFW density profiles with regard to KLAOl simulations. 
This agrees with the results found using standard deviations from the mean, 
represented in Figs. 1-8, which have been deduced directly from simulations, 
with regard to a selected, fitting density profile. Then a valid criterion for the 
choice of a fitting density profile, in connection with a given set of simulated, 
dark matter haloes, appears to be the following. 

Statistic razor for fitting density profiles to simulated, dark matter haloes 

Given two or more fitting density profiles and a set of simulated haloes, the 
best fit is related to the minimum value of the ratio of the absolute error to 
the corresponding standard deviation from the mean, \xrf\/asri = \v~V*\/'^srj- 

4.2 Interpretation in terms of the spherical top-hat 
model 

The spherical, top-hat model makes a valid reference for simulated haloes 
(e.g.. Cole & Lacey 1996; KLAOl), and for this reason we think an interpre- 
tation of the dimensionless parameter, 6, in terms of the spherical top-hat 
model, may be of some utility. With regard to fitting, dark matter haloes, 
the combination of Eqs. (5a), (7), and (8) yields: 

6=^C^E^ P , ^ ■ (31) 

3 ' Miokpc-^Mio ' ^ ' 

which shows the dependence of 5 on the product, pM . 

With regard to FMOl simulations, an inspection of Table 3 shows that 
while the mean density, p200! is inversely proportional to the (fiducial) total 
mass, M200J the dimensionless parameter, 5, is directly proportional to M2oo- 
It can also be seen that the product, P2oo^2oo! also increases as M200 does 
and vice versa, in agreement with Eq. (31). The parameter, 5, appears to 
depend on the mass, M, and a second independent parameter, which may 



38 



be chosen as the mass excess at the start of simulation, {6M/M)start, or the 
related mean peak height, Ustart = i^M/M) start/ < (SM/M) 



start 



> 



1/2 



In the framework of the spherical, top-hat model, the total mass is con- 
served and the following relations hold: 



Pv 



R„ 



Pv 



R^ 



Rn, 



R^ 



(32) 



where the indices, rec, max, vir, denote recombination, maximum expansion, 
virialization, respectively, and the ratio, {Rrec/Rmax)~^, is a solution of the 
third-degree equation: 



1 \.^g(^i/y I S X U Lip 



1 






-A, 



tAj I U Uir' 



1 + 



'6M*' 



; 
(33) 



where ^2 and A represent the density parameter related to matter and cosmo- 
logical constant, respectively, and SM* /M is the mean mass excess within a 
spherical volume where the mass of the material Hubble flow equals M (e.g., 
Lokas & Hoffman 2001a,b). 

The solution to Eq. (33) which is of interest here, in absence of cosmolog- 
ical constant i.e. A ^ 0, has to attain the limiting expression (e.g., Peebles 
1980, Chap.H, §19a): 



Rr 



1 



^^rec 



{5M*/M\ 



{R 



■max) A=0 



red 



\i-rt::A < 



1 + {6M*/M)r 






M 



< 1 



(34) 



where the mass excess, at recombination epoch, has a substantial contribu- 
tion from both the growing and the decreasing mode of the density pertur- 
bation. For further details see e.g., Caimmi et al. (1990). 
The combination of Eqs. (32) and (34) yields: 



Pr. 



{,'n 



'SM*' 



Sn 



{R. 



max) A=0 



|l-^recl < 



'6M*' 



rec 

(35) 
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and, in addition: 



c 



^ — Stiir _ 



(36) 



where the upper hmit corresponds to null kinetic energy at maximum ex- 
pansion, and the lower limit is related to a necessary condition for maximum 
expansion, R^r < Rmax- For further details see e.g., Caimmi et al. (1990). 
Finally, the combination of Eqs. (32), (35), and (36), produces: 



Pvir Smax'^vir\Ph)r 



'6M*' 



11 -fi~!.| < 






<1 



(37) 

where {ph)rec = Prec/[1 + {^M* /M)rec] IS the density of the material Hubble 
flow at recombination epoch. The above relation may be cast under the 
equivalent form: 



Pv 



125 

72^ 
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maxSvtr 
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-\3 



OJ 



|l-^recl< 






<1 



(3J 



where {6M/M)q = (3/5)((5M*/M)rec(l + Zrec) represents the present-day 
mass excess of the growing mode predicted by the top-hat model, in a flat 
cosmology with a vanishing quintessence or, in particular, cosmological con- 
stant. For a formal derivation, see Appendix A. 
The combination of Eqs. (31) and (38) yields: 

3 



5 = CsCn 



M 



C, 



125 



6M' 



y^gi/o'^okpc^' 

o„i — I 



M 



10 



3 



G M 



a,,\\ < 1 



(39a) 



(39b) 
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where the coefficient, Cs, may be evaluated numerically by use of Eqs. (5b), 
(10), and (36), taking the value of the scaled radius, S, from Table 2, and 
keeping in mind that Ho = 50 and 70 km s~^ Mpc~^ have been assumed in 
FMOl and KLAOl simulations, respectively. The result is: 



7.12786 10"^ < iCs)NFW,FM < 5.70229 10"^ 
5.93342 10"^ < iCs)MOA,FM < 4.74673 10"^ 
1.05903 10"^ < iCs)NFW,KLA < 8.47224 10"^ 
8.12553 10~^ < iC5)MOA,KLA < 6.50043 10"^ 



(40a) 
(40b) 
(40c) 
(40d) 
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mass 


FMOl 


KLAOl 


'o^€. 


o-Q 


h 


&2 


log 62 


O-Q 


h 


b2 


log 62 


-1 


15.08 


1.00 10^ 


3.43 10^ 


2.535 


4.200 


1.72 10^ 


7.41 10° 


0.870 





11.85 


2.75 10^ 


1.66 10^ 


3.221 


3.300 


4.73 10^ 


3.59 10^ 


1.555 


1 


8.61 


8.14 10^ 


6.38 10^ 


3.805 


2.400 


1.40 10^ 


1.38 10^ 


2.141 


2 


5.71 


2.64 102 


1.86 10^ 


4.270 


1.590 


4.55 10^ 


4.02 10^ 


2.604 


3 


3.66 


8.89 102 


4.90 10^ 


4.690 


1.020 


1.53 10^ 


1.06 10^ 


3.026 


4 


2.05 


3.42 10^ 


8.61 10^ 


4.935 


0.570 


5.90 10^ 


1.85 10^ 


3.268 


5 


0.937 


1.61 10^ 


8.23 10^ 


4.915 


0.261 


2.77 10^ 


1.77 103 


3.249 


6 


0.345 


9.43 10^ 


4.11 10^ 


4.613 


0.096 


1.62 10^ 


8.84 10^ 


2.947 



Table 11: The mass (in decimal logarithm and unit Mio = IO^^Mq), the 
present-day rms mass excess, ao=< {6M/M)q >^/'^, and the parameters, 
bi = i^recCmaxRniax/^pc, 62 = 6 / {u^^^Cs(rnax) , related to top-hat, sphcrical 
perturbations with same mass spectrum as in Gunn (1987), but normalized to 
cosmological models assumed in FMOl and KLAOl simulations, respectively. 
The plot of the mass spectrum (Gunn 1987) has been assumed to reproduce 
only the growing mode. 



in connection with both NEW and MOA density profiles, related to both 
FMOl and KLAOl simulations. 

With regard to FMOl (as = 0.7) and KLAOl (as = 0.9) simulations, the 
mass enclosed within a spherical region of radius Rg = 8h~^ Mpc, is: 



{Ms)fm = 

iMs)KLA 



2.021868 lO^Mio 
= 4.842291 lO^M 



10 



(41a) 
(41b) 



for a formal derivation, see Appendix A. The normalization of top-hat, spheri- 
cal perturbations with same spectrum as in Gunn (1987), to the above values, 
demands a multiplication of the rms mass excess plotted in Gunn (1987) by 
a factor of 14/13 and 3/10, respectively. The related values, together with 
some other parameters which are characteristic of top-hat, spherical pertur- 
bations, are listed in Table 11. The parameter, S, expressed by Eqs. (39), is 
plotted as a function of the mass, for lower and upper values of the factor, 
Cs, expressed by Eq. (40), and mean peak heights Urec = 1,2,3,4, in Fig. 9, 
where the values deduced from FMOl (top panel) and KLAOl (bottom panel) 
simulations, in connection with MOA (top panel) and NFW (bottom panel) 
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density profiles, respectively, are also represented. In both cases it is ap- 
parent that, for increasing masses, simulated haloes correspond to fitting, 
top-hat haloes with increasing (mean) peak height. 

The occurrence of cosmological constant changes the value of the treshold, 
beyond which density perturbations are expanding forever (e.g., Lokas & 
Hoffman 2001a,b). The mass of the critical density perturbation, related to 
an infinite turnaround radius, is marked by a vertical line in Fig. 9, from the 
left to the right for increasing peak height. 

Keeping in mind that, for a fixed peak, upper curves are related to no 
acquisition of angular momentum {(mr = 8) and lower curves are related 
to maximal, allowed acquisition of angular momentum {(mr = 1), the corre- 
spondence between simulated and fitting, spherical, top-hat haloes implies an 
efficient acquisition of angular momentum, especially for low masses. On the 
other hand, with regard to FMOl runs, simulated haloes - especially at low 
masses - appear to be consistent with fitting, spherical, top-hat haloes with 
mean peak heights within the range 1 ~ Urec ~ 2, contrary to what would 
be expected, Urec ~ 2; with regard to KLAOl runs, the above mentioned 
discrepancy disappears. 

5 Concluding remarks 

Simulated dark matter haloes have been fitted by self-similar, universal den- 
sity profiles, where the scaled parameters depend only on a scaled (trun- 
cation) radius, S = R/tq, which, in turn, has been supposed to be inde- 
pendent of the mass and the formation redshift. The further assumption of 
a lognormal distribution (for a selected mass bin) of the scaled radius, or 
concentration, in agreement with the data from a large statistical sample of 
simulated haloes (Bullock et al., 2001), has allowed (at least to a first extent) 
a normal or lognormal distribution for other scaled parameters, via the same 
procedure which leads to the propagation of the errors. 

A criterion for the choice of the best fitting density profile has been pro- 
posed, with regard to a set of high-resolution simulations, where some aver- 
aging procedure on scaled density profiles has been performed, in connection 
with a number of fitting density profiles. To this aim, a minimum value of the 
ratio, \xrj\/asrj = \v ~ V*\/'^srj, has been required to yield the best fit, where 
r] is the arithmetic mean over the whole set; r]* is its counterpart related to 
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Figure 9: The dimensionless parameter, S, vs. the fitting mass, M, for 
top-hat, spherical perturbations, with no (open circles) and maximum al- 
lowed (crosses) acquisition of angular momentum, according to Eqs. (40b) 
and (40c), in connection with MOA (top panel) and NFW (bottom panel) 
density profiles. For each mass and class of symbols, values from down to 
up are related to mean peak heights at recombination epoch, Urec = 1, 2, 3, 4, 
respectively. The mass of the critical density perturbation, for which the 
turnaround radius is infinite, is marked by a vertical line, from the left to 
the right for increasing peak height (bottom panel). Values related to FMOl 
(top panel) and KLAOl (bottom panel) simulations are also represented, 
with same symbol captions as in Figs. 1-4 and 5-8, respectively. Passing from 
NFW to MOA density profile for each set of simulations, would lower values 
related to top-hat, spherical perturbations, by about one dex, and vice versa. 
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the fitting density profile; asjj is the standard deviation from the mean; and 
?7 is a selected, scaled i.e. dimensionless parameter. 

The above criterion has been applied to a pair of sets each made of a dozen 
of high- resolution simulations, FMOl (Fukushige & Makino 2001) and KLAOl 
(Klypin et al. 2001), in connection with two currently used density profiles, 
NFW (e.g., Navarro et al. 1997) and MOA (e.g., Moore et al. 1999), where 
the dependence of the scaled radius on the mass and the formation redshift, 
may be neglected to a first extent. With regard to FMOl and KLAOl samples, 
the best fit has revealed to be MOA and NFW, respectively. In addition, the 
above results have been found to hold also in dealing with rms errors derived 
via the propagation of the errors, with regard to the distributions of scaled 
parameters. The sensitivity error of simulations has also been estimated and 
shown to be less than the related, standard deviation, that is a necessary 
condition for detectability of accidental errors. 

Some features of the early evolution of dark matter haloes represented 
by fitting density profiles, have been discussed in the limit of the spherical 
top-hat model. Though the related matter distributions have appeared to be 
poorly representative of simulated haloes, unless the (mean) peak height is an 
increasing function of the mass, the results have been shown to be consistent, 
provided considerable acquisition of angular momentum took place during 
the expansion phase. 
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7 Appendix 

7.1 A. Some properties of spherical, top- hat density 
perturbations 

With regard to a special class of QCDM cosmological models, where quintessence 
obeys the equation of state p* = wp^l relating its density, p^, to its pres- 
sure, p.1,, via a time-independent parameter or quintessence index, w, where 
— 1 < w < 0, the following relations hold (e.g., Lokas 2002; Lokas & Hoffman 
2002): 
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where a is the scale factor normalized to unity at present, z is the redshift, fi 
and \E' the density parameter of matter and quintessence, respectively, H the 
Hubble parameter, and the index denotes the current time. In the limit of a 
vanishing quintessence, Eqs. (42), (44), and (46) reduce to their counterparts 
in cosmological models with sole matter and radiation (e.g., Zeldovich & 
Novikov 1982, Chap. HI, §4). The special case, w = —1, corresponds to the 
cosmological constant. 

An equivalent expression of Eqs. (44) and (45), using (46), is: 
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and Eq. (47) allows the validity of the inequality: 

' ' "^ <Q-Uz)<l + - — ^ ; -l<w<—. 



1 + z Qo + '^o ~ 1 + 2 f^o --3' 

(49) 

which shows that the evolution of the matter density parameter in a cos- 
mological model defined by assigned values of (rioj^o)) is comprised be- 
tween its counterparts related to cosmological models defined by (Qq, 0) and 
{Qq + \E'o, 0), respectively. 

With regard to the function, 0(a), defined by Eq. (43), a detailed analysis 
involving the solution of a third-degree equation yields the final result: 



1+^0 



^^^ < Ma) < 



1 I ^0 /^ -1 1 

1 H [a — 1 



-1/2 



Qn > 0.25 



2 

(50) 

which shows that the evolution of the scale factor in a cosmological model 
defined by assigned values of {Qq, "^q), is comprised between its counterparts 
related to cosmological models defined by {fio/2, 0) and {Qq, 0), respectively, 
provided Qq > 0.25. 

With regard to the Hubble parameter, H(t), defined by Eq. (46), a de- 
tailed analysis involving the solution of a transcendental equation yields the 
final result: 

hl + z)\l + noz)<^^-^<{l + z)\l + noz) ; f]o>0.25 ; (51) 

which shows that the evolution of the Hubble parameter in a cosmologi- 
cal model defined by assigned values of (fig, \E'o, -f^o), is comprised between 
its counterparts related to cosmological models defined by (f2o, 0, Hq) and 
(fio, 0, ii^o/v^), respectively, provided VLq > 0.25. 

Cosmological models with sole matter and radiation evolve at the same 
rate in the limit |1 — Q~^\ <^ 1 i.e. at early times (e.g., Zeldovich & Novikov 
1982, Chap. HI, §4) and, owing to inequalities (49), (50), and (51), the same 
holds for cosmological models with quintessence. In the above mentioned 
limit, the matter density of the Hubble flow, ph, reads (e.g., Peebles 1993, 
Chap. II, §13): 

where G is the constant of gravitation. 
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The combination of Eqs. (37) and (52) yields: 
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<1 



(53) 



in the framework of the top-hat, spherical model, related to a flat universe 
with vanishing quintessence, the growing mode of the density perturbation 
attains the present-day value (e.g., Peebles 1980, Chap. II, §15): 

^6M\ 3 f6M*\ ,^ 

■ (54) 
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and the combination of Eqs. (53) and (54) yields Eq. (38). 

In addition, the combination of Eqs. (34), (35), (52), and (54) produces: 
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(55) 

where the mass excess has been expressed as the product of the present-day 
rms mass excess and the peak height at recombination epoch, averaged over 
the whole volume, that is: 
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for more details see e.g., Caimmi et al. (1990). Of course, Eq. (55) makes 
a lower limit to the product, i^recRmax, as Eq. (34) holds in the limit of a 
vanishing quintessence. 

On the other hand, in open or flat universes, density perturbations below 
a treshold are destined to expand forever. In the special case of cosmological 
constant, the critical value is (e.g., Lokas & Hoffman 2001a,b): 
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where the turnaround occurs at an infinite radius. 

With regard to FMOl and KLAOl simulations, the present-day rms mass 
excess in a spherical region of radius Rs = 8h^^ Mpc, takes the value: 
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owing to the general definition of mass excess: 

SM p-ph 



M 



Ph 



(59) 



the mass within the region under consideration may be obtained using Eqs. (52), 
(58), and (59). The result is: 



Ms = 2^ 10^ (1 + as) 






(60) 



which may be specified for any fiat cosmological model with assigned values 
of iJo, ^0, and ag. 

7.2 B. Determination of fitting haloes to FMOl and 
KLAOl simulated haloes, with regard to NFW and 
MOA density profiles 

Given a selected, FMOl simulated halo among the runs listed in Table 3 i.e. 
for which the values of the parameters, M2005 ''"200 > and 5, are known, we 
aim to derive the values of the scaling density, po, and the scaling radius, 
ro, related to the corresponding, fitting halo, and then the remaining param- 
eters, in connection with both NFW and MOA density profiles. The total 
mass of the simulated halo, M, appears in FMOl prescriptions, expressed by 
Eqs. (4) and (5), but the related values are not reported therein. Then the 
key parameter is the scaled radius, ^200- 
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Owing to the general definition of scaled radius, Eq. (5a) may be written 
under the equivalent form: 



e 



200 



^200 
^0 



kpc 



M 



M 



10 



-2/3 



(61) 



on the other hand, the particularization of the general expression of the mass 
enclosed within a generic, isopycnic surface, to the case under discussion, via 
Eqs. (2) and (3), reads: 



M2oo = 3Mo/^'°V(Oe'de 
Jo 



(62) 



and the combination of Eqs. (61) and (62) produces a transcendental equation 
in M, which can be solved in connection with an assigned (NEW or MOA) 
density profile, provided the values of the parameters, C^, 5, r2oo, and M200, 
are specified using the data listed in Table 3. 

The knowledge of the total mass of the simulated halo, assumed to co- 
incide with the mass of the fitting halo, M, allows the calculation of the 
scaling density, po? and the scaling radius, tq, via Eqs. (4) and (5), and then 
the radius along a fixed direction, R = Srg, the scaling mass, Mq, appearing 
in Eq. (6), and the dimensionless parameter, k, expressed by Eq. (8). 

With regard to NEW density profiles, the particularization of Eq. (62) to 
the case under discussion, yields: 



M 



200 



12Mn 



ln(l + e 



e 



200 



200 J 



,200 



(63) 



where, owing to Eqs. (6) and (19): 

^Maoo ^ 125 M200 
12 Mo 287r M ' 

and the combination of Eqs. (63) and (64) yields: 



(64) 



M 



200 



M 



287r 
125 



e 



In 



1 



200 



1+6 



200 



(65) 



finally, the combination of Eqs. (5b), (61), and (65) produces the ultimate 
transcendental equation in M. 
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With regard to MOA density profiles, the particularization of Eq. (62) to 
the case under discussion, yields: 

M2oo = 4Mo In (1 + ^20^0') ; (66) 

where, owing to Eqs. (6) and (17): 



1 M200 375 M200 



4 Mo 567r M ' 
and the combination of Eqs. (66) and (67) yields: 



(67) 



200 ODTT / 3/2\ .„„v 

— = In 1 + ^200 ; (68) 



^^200 _ 567r 
lir"375^"V^^''2ooy 

finally, the combination of Eqs. (5b), (61), and (68), produces the ultimate 
transcendental equation in M . 

The above procedure, via Eqs. (10)-(16), also holds for a selected, KLAOl 
simulated halo among the runs listed in Table 4, for which the values of the 
parameters, Mmr-, f^vir, and S, are known. 

With regard to both NEW and MOA density profiles, some parameters 
related to fitting haloes, in connection with simulated haloes from FMOl and 
KLAOl, are listed in Tables 7 and 8, respectively. 

7.3 C. Sensitivity errors of dark matter halo simula- 
tions 

Bearing in mind the general results listed in Tab. 1, together with Eq. (8), 
the scaled parameters, M/Mq, R/tq, p/po, and k, depend only on the fitting 
density profile. For a selected choice of exponents (q;,/9,7), fitting haloes 
depend on two parameters, (ro,po)) or {M,6), via Eqs. (4) and (5). 

Given a computer run with A^ identical particles of mass m, the sensitivity 
error with respect to the mass is clearly expressed as AM = m. It follows 
that: 

A ^200 A M2oo \ M200 m 

^Mr=l^ + l^j mTa^ ' ^^o = AM2oo = m ; (69) 
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the second parameter, 6, is proportional to the present-day mass excess of 
the growing mode predicted by the top-hat model, {SM/M)o, as: 



S = C, 



'SAT 
M 



M 
Wo ' 

2o 1.^^3 



C, 



125 g^a i^o^l^okpc 

54^"^ '•" G Mio 



(70a) 
(70b) 



where Qq is the present-day, matter density parameter (r2o + Aq = 1), and Cmr 
depends on the evolution of the density perturbation during the expansion 
phase, and lies within the range 1 < (mr < 8. For a formal demonstration, 
see Appendix A. 

The repetition of the above procedure yields: 



.^ 3M2oo + 2(5M)o ^ 

^^ = ^ — 7-? — 7Tl^-^^ "^ = <J 

M2oo(5M)o 

A(5M)o = AM200 = m ; 



2 + 3Cy'6-'/' (^ 
' V Mio 



1/3' 



m 



M 



200 



;(71a) 
(71b) 



and the sensitivity error with respect to po and Tq, by use of Eqs. (4) and (6), 
is: 



M200 + (^M)o 
Apo = Spo^TT 7777^ "^ = 3P0 



1 
Aro = 0^0 



M2oo(5M)o 

4 + 3cy^(5-i/2 



1 + cl/3^-1/3 (^ 



M 
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M 



10 
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i(72) 
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finally, the further assumption Ar 
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M200 
Atq allows the following results: 
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(74) 
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where Eqs. (4) and (7) have been used. 

The sensitivity error of FMOl computer runs in dark matter halo simu- 
lations, expressed by Eqs. (69), (74), (75) and (76), may be calculated and 
compared with their rms counterparts, to see if a necessary condition for the 
detectability of accidental errors is satisfied, namely: 

Ar]<(r^ ; (77) 

where 7] = M200/M0, r2oo/^o, P200/P0, /«200- To this aim, an inspection of 
Tabs. 3, 5, and 7, shows that the following inequalities hold: 

—— < lO"*' ; — ^ < 20 ; ^^ < 20 

M200 Mo ' ro 

^<1 ; .200 < 2.8 ; C^VV^ f^V' 
Po V Mio 

where upper values of Cs have been used, in connection with the range 1 < 
Cvir < 8, according to Eq. (39b). The combination of Eqs. (69), (74), (75), 
(76), and (78) yields: 

Mo ro Po 

(79) 

The sensitivity errors of KLAOl computer runs in dark matter halo sim- 
ulations, are expressed in the same way, provided ^200 is replaced by u^ir 
therein, where u = M, r, p, k. An inspection of Tabs. 4, 6, 8, shows that the 
following inequalities hold: 

^ < 6 10-^ ; ^ < 30 ; ^ < 20 ; ^ < 0.2 ; (80) 



^200 2 


(78a) 
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<0.2 ; 


(78b) 



M^ir Mo ro r. 
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Pvir ^ ^o ^ on /^l/^r— 1/3 ( ^^^vir\ 



<0.02 ; /€,„;.< 20 ; cT^d-^'H^] < 1.2 

10 



Po VM 



where upper values of Cs have been used, in connection with the range 1 < 
Cvir < 8, according to Eq. (39b). Following the same procedure used for 
FMOl, yields the final result: 

A^<1 ; A^<4 10-^ ; A^<5 10-^ ; A/s;,,, < 6 10^^ . 



Mo ro Po 

(82) 
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The comparison between the sensitivity errors, expressed by Eqs. (79) 
and (82), and their rms counterparts, deduced from Table 9, shows that a 
necessary condition for the detectabihty of accidental errors, expressed by 
Eq. (77), is satisfied for both FMOl and KLAOl simulations. To this aim, it 
is worth remembering that o-j^ = (T^/\/iV, where iV = 12 in the case under 
consideration, according to the theory of the errors. 

7.4 D. rms errors of distributions depending on scaled 
parameters 

Let dark matter haloes be fitted by universal density profiles, expressed by 
Eq. (2), and let the distribution depending on the scaled radius, S (or con- 
centration with regard to NEW density profiles), be lognormal. The scaled 
mass enclosed within the generic scaled distance, ^, and the related scaled 
mean density, are: 



M{0 



Mr 







P(0 3 f^ 



3^ /(Oe'de ; (83) 

Jo 

tfiOed^ ; (84) 

Jo 



Po Uo 

and the generalization of the dimensionless parameter, k, defined by Eq. (8), 
to the generic scaled radius, ^, reads: 

m = ci'^e" ; (85) 

where the constant, C^, is determined by averaging on the results of simula- 
tions (EMOl), and for the cases of interest it is expressed by Eqs. (5b) and 
(10). 

The first derivatives of the functions on the left-hand side of Eqs. (83), 
(84), and (85), are: 

d(M/Mo) .,,_, . ^3^^ 

(87) 
(88) 
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and, in addition: 

dioge 1 1 



(89) 



de InlOe 

following the standard rules of derivation. 

Let us suppose that (i) the scaled parameters, rj, rj = Mtrn/Mo, Etrn, 
Ptrn/po, and ktrn, tm = 200, vir, as functions of logEtm, are expressible as 
Taylor series where the starting point coincides with the expected value of 
the lognormal distribution of the concentration, S = exp^Q(\ogEtrn)*] (h) the 
convergence radius of the series under discussion exceeds at least three times 
the rms error of the lognormal distribution, (JiogHt™) or in other words the 
convergence occurs at least within the range, logS=F3(TiogHt^„; (iii) the series 
under discussion can safely be approximated by neglecting all the terms of 
higher order with respect to the first. It is worth remembering that the prop- 
agation of the errors, and the related formulae currently used in literature, 
are grounded on the above mentioned assumptions. 

Let 0(^) be a generic, derivable function of an independent variable, ^. 
The application of the theorem of the derivative of a function of a function, 
where the second function is log^, yields: 

dlog^ d^ dlog^ d^ 

and the particularization to the scaled parameters considered here, reads: 



dlog^ 

d(p/po) 
dlog^ 



3 In 10 /(Or ; (91) 



3 In 10 



^^^' e^ Mo 



(92) 



:^ In 10 MO ; (93) 



dloge 2 

owing to Eqs. (86), (87), and (88). 

On the other hand, the validity of the above assumptions implies that 
the relations: 

rj{\ogEtrn) = ^(logS) + — — - — (logSt^n - logS) ; (94a) 

; trn = 200, vir ; (94b) 



Mo' ^' 7o'' 
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hold to a good extent. The combination of Eqs. (89), (91), (92), (93), and 
(94), yields: 



M(log : 
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Po 

3 
A; + -lnlO(logE:t^„-logE:) ; 



(logEtrn -logS) 



(95) 
(96) 
(97) 

(98) 
(99) 



where M = M(S), p = p(S), and k = A;(S), for sake of brevity. 

The scaled parameter, logSf^n, niay be considered as a physical quantity 
to be measured directly. Accordingly, the distribution depending on logSj^n 
has necessarily to be normal. A theorem related to the theory of errors 
ensures that the scaled parameters, 77, defined by Eq. (94), also follow normal 
distributions, whose expected values and rms errors, via Eqs. (95), (96), (97), 
and (98), are expressed as: 
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O'logH 



(100) 
(101) 



and the last relation, owing to Eqs. (89), (91), (92), and (93), takes the 
explicit form: 



o-Af/Mo = 3 In 10 S'V(2) o-iogH ; 
(Js = InlO H (TiogH ; 

1 M 



(j-plp^ =3 In 10 /(H) 
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23 Mn 
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lg= ) 



(102) 
(103) 

(104) 
(105) 



for the scaled parameters under consideration. 

Starting from Eq. (101), after replacing r] with log?], and using Eq. (89), 
after replacing ^ with r], yields an expression of the rms errors of lognormal 
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distributions, in terms of their counterparts related to normal distributions. 
The result is: 
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or, after comparison with Eq. (101): 
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which, owing to Eqs. (102), (103), (104), and (105), take the explicit form: 
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where the identity has been written for sake of completeness. 
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7.5 E. Random model: the concentration distribution 
as a result of the central limit theorem 

Dark matter halo and star formation take place in a similar fashion, namely 
a transition from an undifferentiated fluid to substructures. Though a molec- 
ular cloud is neither expanding nor subjected to the Copernican principle, 
contrary to the Hubble flow, still the above mentioned processes are expected 
to exhibit some common features. 

The initial mass function in a star generation may safely be fitted by a 
lognormal distribution (e.g., Adams & Fatuzzo 1996; Padoan et al. 1997), 
which, in turn, can be interpreted in terms of the central limit theorem 
(Adams & Fatuzzo 1996). On the other hand, data from a statistical sample 
of about five thousands of simulated dark matter haloes (Bullock et al. 2001), 
show - to a good extent - a lognormal distribution of the concentration within 
mass bins of (0.5-1.0) xIO^/i^-'^Mq, where 11 < n < 14 and n is an integer. 



Here we adopt a statistical approach to the calculation of the lognor- 
mal distribution of the concentration, following the same procedure used by 
Adams & Fatuzzo (1996), in dealing with the initial mass function of a star 
generation. To this aim, first let us suppose that a transformation exists 
between initial conditions and the final properties of the dark matter halo 
with assigned mass. 

Given a cosmological model and a perturbation spectrum, the initial con- 
ditions of a simulation are defined by a generation of complex numbers with 
a phase randomly distributed in the range < < 27r and with ampli- 
tude normally distributed, where the variance is provided by the selected 
spectrum, in the simplest case of a Gaussian distributed random field. For 
further details see e.g., Moscardini (1993); Tormen et al. (1997). 

Given a universal density profile, assumed to fit to dark matter haloes 
under consideration, the final properties are related to the scaled radius, or 
the concentration with regard to NFW density profiles, via the results listed 
in Table 1. 

Second, let us suppose that the transformation under consideration is 
expressible as a product: 



■^trn 



^n/5? = n«. ; (112) 



3 

i=i i=i 



where the constant. A, and the exponents, 7^, are fixed, and the variables, 
l3j or aj, are conceived as random variables. Though Eq. (112) cannot be 
motivated by the existence of a semiempirical relation, as in the case of star 
formation (Adams & Fatuzzo 1996), neverthless it cannot be excluded unless 
further knowledge about the genesis of dark matter haloes will be available. 

The central limit theorem holds provided the random variables, aj, ap- 
pearing in Eq. (112), are completely independent, and their number, n, tends 
to infinite. For the more realistic case of a finite number of not completely 
independent random variables, the resulting distribution is expected to be 
different from a (log) normal distribution. 

Taking the decimal logarithm on both sides of Eq. (112) yields: 

n 

log St„, = ^ log aj ; (113) 

and the expected value of the distribution, fj(\ogaj) dloga^, depending on 
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the random variable, log aj, is: 

/ + 00 
^ogajfj{logaj)dlogaj ; I < j < n ; (114) 

-oo 

accordingly, the error, Xioga , is: 

xioga, = logctj - (logftj)* = log^ ; (115a) 

aj 

a* = expio(logaj)* ; 1 < j < n ; (115b) 

where, of course, a* is different from the expected value of the distribution, 
fj{oij) dttj, depending on the random variable, aj. 

The variance of the distribution, fj{x\ogaj) dxiogo^, is: 

/+00 
a^foga,/i(a;ioga,)dxioga, ; l<j<n ; (116) 

-oo 

and the related, expected value, equals zero. 
Let us define the random variable: 

n n n n 

C = E ^log". = E log -i = E log "i - E log "i 5 (117) 

and combine Eqs. (113) and (117), to obtain: 



iogSt,„ = c + Eiog"j* ; (118) 



which is equivalent to: 



Str-n = S*expio(C) ; (119a) 

n 

s* = n«; ; (119b) 

where, of course, S* is different from the expected value of the distribution, 
f(Etrn) dSirn, depending on the random variable, St™. 

The application of the central limit theorem to the distribution, f{Q d(, 
depending on the random variable, (, yields (e.g., Adams & Fatuzzo 1996): 

n 

^? = E<a. ; (120) 
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and aiming to find a distribution, f{() d(, characterized by unit variance, 
a^ = 1, and null expected value, (* = 0, let us define the random variable: 

C = 4 ; (121) 

where, owing to the central limit theorem, the distribution is normal: 

fiC) dC = -^exp(-^)dC ■ (122) 



independent of the initial distributions, /j (xiog « J dxiog q^. . 

Taking the decimal logarithm of both sides of Eq. (119a), and using 
Eq. (121), yields: 

log Etrn = log S* + C = log E* + ale ; (123) 

and the distribution, fiEtm) dEtm, depending on the random variable, Etm, 
is lognormal. The related expected value and variance are: 



(logSi„J* 


= logS* ; 


(124) 


<s.r. = ^! ; (125) 


accordingly, the distribution reads: 




r (logS,,„-logS*)21 
2^2 


dlogSi,„(^26) 


J [log ^trn ) 0. log ^trn — /-— CXp 
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and the decimal logarithm of the probability density, /(logS^^^), may be 
written as: 



i , /„ N , i i / , I — 'trn 



log[/(logS,,„)] = -^log(27r)-logac-j^^(^log^j' ; (127) 

where the first term on the right-hand side member may be conceived as a 
normalization constant. 

The values of the expected value and the variance, expressed by Eqs. (124) 
and (125), related to the lognormal distribution, defined by Eq. (126), may 
be deduced from the results of simulations (e.g., Bullock et al. 2001). For 
further details on the procedure outlined above, see Adams & Fatuzzo (1996). 
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By analogy with the theory of measure, a computer run may be consid- 
ered as an execution of measure operations, the related computer code as a 
measure instrument, the dark halo as a measure subject, and the sequences 
of random numbers used in the definition of initial conditions as contributors 
to the accidental error. Then the computer output may be thought about as 
a measure of the corresponding scaled parameter, which may be conceived 
as fluctuating around its fitting counterpart. 

In addition, it is worth of note that the application of a least-squares 
or least-distances method in fitting simulated with universal density profiles 
(e.g., Dubinski & Carlberg 1991; KLAOl; FM03) implies a (fiducial) normal 
distribution of the simulated density (in decimal logarithm) around the ex- 
pected value deduced from the fitting density profile, at any fixed distance 
(in decimal logarithm). It is the particularization, to the case of interest, of 
a well known result of the theory of errors (e.g., Taylor 2000, Chap. 8, § 8.2). 
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